arXiv:nucl-th/9706060vl 24 Jun 1997 


Fusion barrier distributions 
in systems with finite excitation energy 


K. Hagino , 1,3 N. Takigawa , 1,3 and A.B. Balantekin 2,3 


1 Department of Physics, Tohoku University, Sendai 980-77, Japan 

2 Physics Department, University of Wisconsin, 

Madison, Wisconsin 53706 

3 Institute for Nuclear Theory, University of Washington, 
Seattle, Washington 98915 


Abstract 

Eigen-channel approach to heavy-ion fusion reactions is exact only when the excitation 
energy of the intrinsic motion is zero. In order to take into account effects of finite 
excitation energy, we introduce an energy dependence to weight factors in the eigen- 
channel approximation. Using two channel problem, we show that the weight factors are 
slowly changing functions of incident energy. This suggests that the concept of the fusion 
barrier distribution still holds to a good approximation even when the excitation energy of 
the intrinsic motion is finite. A transition to the adiabatic tunneling, where the coupling 
leads to a static potential renormalization, is also discussed. 
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It has been well recognized that cross sections of heavy-ion fusion reactions at energies 
near and below the Coulomb barrier are strongly influenced by coupling of the relative 
motion of the colliding nuclei to nuclear intrinsic motions|l|. When the intrinsic degree of 
freedom has a degenerate spectrum, fusion cross sections can be calculated in the sudden 
approximation. In this limit, the coupling gives rise to a distribution of fusion potentials. 
The fusion cross section is then given by an average over the contributions from each 
fusion barrier with appropriate weight factors^, [f f|]. Based on this idea, a method was 
proposed to extract the barrier distribution directly from the fusion excitation function 
using the second derivative with respect to the product of the fusion cross section and the 
center of mass energy, Ea [Q. This stimulated precise measurements of the fusion cross 
sections for several systems^]. The extracted barrier distributions have then been found 
to be very sensitive to the structure of the colliding nuclei, e.g. the sign of hexadecapole 
deformation parameter, while the fusion cross section itself is rather featureless|| |7j. 

In a rigorous theoretical interpretation, the barrier distribution representation, i.e. 
the second derivative of Ea , has a clear physical meaning only if the excitation energy 
of the intrinsic motion is zero. Nonetheless, this analysis has been successfully applied 
to systems with relatively large excitation energies || |], |J. For example, the second 
derivative of Ea for 16 0 + 144 Sm fusion reaction clearly shows the effects of coupling 
to the octupole phonon state in 144 Sm, whose excitation energy is 1.8 MeV, much more 
clearly than the fusion cross section itself||. |j. Also the analysis of the fusion reaction 
between &8 Ni and 60 Ni, where the excitation energies of quadrupole phonon states are 1.45 
and 1.33 MeV, respectively, shows that the barrier distribution representation depends 
strongly on the number of phonon states included in coupled-channels calculations ||. 
These analyses suggest that the representation of fusion process in terms of the second 
derivative of Ea is a very powerful method to study the details of the effects of nuclear 
structure, irrespective of the excitation energy of the intrinsic motion. 

Despite these successes, there remains the question whether the second derivative of 
Ea represents a “distribution ”of fusion potential barriers when the excitation energy of 
the intrinsic motion is not zero. To address this problem, one may attempt to use the con¬ 
stant coupling approximation, where the coupling form factor is assumed to be constant 

m 


throughout the interaction range 


In this approximation, the unitary transformation 


gives the fusion cross section as a sum of eigen-channel cross sections. The constant cou¬ 
pling approximation has, however, a serious drawback in that the results strongly depend 
on the position where the strength of the coupling Hamiltonian is evaluated. Therefore, 
one needs to use a different approach to address the question. 

In this paper we treat the weight factors as energy dependent variables. The possi¬ 
bility of the energy dependence of weight factors when the excitation energy is finite has 
been suggested in ref. 0- Here, we explicitly study the energy dependence by perform¬ 
ing exact coupled-channels calculations. It will be shown that the energy dependence 
is quite weak for a wide range of excitation energy, suggesting that the eigen-channel 
approximation works well even when the excitation energy of the intrinsic motion is not 
small. 

Let us consider, for example, the case where the intrinsic degree of freedom has only 
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two levels. The coupled-channels equations then read 
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where U(R) and /i are the bare potential and the reduced mass of the relative motion, 
respectively. F(R) is the coupling form factor and e is the excitation energy of the 
intrinsic motion, respectively. When the excitation energy e is zero, the unitary matrix 
which diagonalizes the coupling matrix is independent of the position R , and the barrier 
penetrability is exactly given by 

P{E) = l - [Pq(E; U(R ) + F(R)) + P 0 (E; U{R ) - F(R))] , (2) 


where Pq(E ; V(R)) is the penetrability of potential V(R) at energy E. In this limit, the 
original single barrier splits into two effective potential barrier with equal weights. 

If the excitation energy e is not zero the unitary matrix depends explicitly on the posi¬ 
tion R. The unitary matrix, therefore, does not commute with the kinetic energy operator, 
and hence simple eigne-channel approach does not hold. An approximation which is often 
made is to assume that the coupling form factor F(R) is a constant. In this approxima¬ 
tion, the unitary matrix becomes independent of R, and the penetrability is given by a 
formula similar to eq. (2) but with different eigen-values and weight factors flOj . In order 
to take into account the radial dependence of the coupling form factor, an approximation 
was proposed, where the coupling matrix is diagonalized to obtain the eigen-barriers at 
each position of the relative motion R while the weight factors are calculated at a chosen 
position R n; j [T3| |. The computer code CCMOD uses this approximation |[TJ . 


Even when we introduce this approximation, the weight factors are still functions of 
the chosen position R, w , and the results might strongly depend on that choice. Usually R w 
is chosen to be the position of maximum of the bare potential barrier |T3|, Q. However, 
there is no theoretical justification that this is the optimum choice, and furthermore, it 
is not obvious whether one can determine the weight factors independent of the incident 
energy. 

In order to avoid these drawbacks and examine the energy dependence of the weight 
factors, we parametrize the penetrability as 


P(E) = v + (E)P + (E) + v.(E)P_(E), (3) 

and evaluate the weight factors v± at each incident energy E. Here P±(E) are the pene¬ 
trabilities of the eigen-potentials U(R) + \±(R), respectively, A ±(R) being the eigen-values 
of the coupling matrix at each position of R , which are given by 

A ±(R) = (e ± \] e 2 + 4 F{R)^j /2. (4) 

If the weight factors slowly vary as functions of the incident energy, the effects of channel 
coupling can be interpreted in terms of the barrier distribution, even when the excitation 
energy is non-zero. In the appendix, we generalize eq. (3) to multi-channel case using the 
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path integral approach. Since the weight factors have to satisfy the nnitarity condition^], 
they can be uniquely determined in the present two level problem, and are given by 


v + (E) = (P(E)-P4E))/(P + (E)-P_(E)), (E) = (P + (E)-P(E))/(P + (E)-P_(E)). 

(5) 

We performed coupled-channels calculations to examine the energy dependence of the 
weight factors. To this end, the coupled-channels equations have to be solved with good 
accuracy, since the penetrabilities P(E) and P±(E) in eq. (5) are exponentially small 
quantities at energies below the barrier. The incoming wave boundary condition, which 
is often used in coupled channels calculations for heavy-ion collisions |H| , can bring some 
numerical errors, though they would be small enough for the purpose of calculating fusion 
cross sections. In order to avoid this, we use here a schematic model of heavy-ion fusion 
reactions by Dasso et al. |1(J, where the radial motion between colliding nuclei is treated 
as a one dimensional barrier penetration problem. Following ref. [10 , we assume the 


gaussian shape for both the bare potential and the coupling form factor, 

U(R) = U 0 e- R2/2s \ F(R) = F 0 e~ R2/2s2 f. 


( 6 ) 


The parameters are chosen following ref. [|l(J to be f/ 0 = 100 MeV, F 0 = 3 MeV, and 


s = Sf = 3 fm, respectively, which mimic the fusion reaction between two 58 Ni nuclei. We 
have checked that our conclusions do not significantly change as long as the value of sj is 
not too small. The mass /i and the excitation energy e are taken to be 29 m n, rnjy being 
the nucleon mass, and 2 MeV, respectively. 

The upper panel of Fig. 1 shows potential barriers for the present problem. The dotted 
line is the bare potential U ( R ) and the solid and the dashed lines are the eigen-potentials 
U(R) + X-(R) and U(R) + A +(R), respectively. If we adopt the eigen-channel picture, 
Fig. 1 means that the potential barriers are “distributed ’’and the original bare potential 
(the dotted line) splits into the two potentials (the solid and the dashed lines). Because 
of the non-commutatibity of the unitary matrix which diagonalizes the coupling matrix, 
and the kinetic energy operator, the standard eigen-channel picture does not apply. If 
we ignore this non-commutatibity (i.e., if we adopt with the CCMOD prescription) the 
weight factors are given by 


w ± (R) = F(R) 2 / (. F(R) 2 + A ± (i?) 2 ) . (7) 

They are shown in the lower panel of Fig. 1 as a function of the distance R. Although they 
do not vary so much near the barrier position, their changes are appreciable throughout 
the barrier region. Therefore, the prescription to fix the weight factors to the values at 
the barrier position may not be satisfactory. 

The results of exact coupled-channels calculations are shown in Fig. 2. The first 
and the second panels are the penetrabilities and their first derivative with respect to 
the energy, respectively. The latter corresponds to the second derivative of Ea ||. We 
use the point difference formula with A E = 2 MeV to obtain the first derivative of 
the penetrability, as is often done in the analyses of heavy-ion fusion reactions^]. The 
first derivative of the penetrability has a clear double-peaked structure, which could be 
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associated with the two eigen-potential barriers and could thus be interpreted in terms 
of a ‘barrier distribution In order to see whether this is the case, we plot in the last 
panel of Fig. 2 the optimum weight factors defined by eq. (5). The solid and the 
dashed lines in the figure correspond to the weight factors for the lower and the higher 
potentials, respectively. We observe that the optimum weight factors change only slightly 
as functions of the incident energy; their change is barely 2.6 % from 10 MeV below the 
barrier to 10 MeV above the barrier. Another important result of this calculation is that 
the weight factors are considerably different from those estimated in the usual way, i.e. at 
the barrier position. When the constant coupling approximation was first introduced, it 
was expected that the choice of the position where the weight factors are estimated is not 
critical to calculate tunneling probabilities because the weight factors vary slowly across 
the barrier region |T0|, |T3| . The position was then chosen to be the barrier position of the 
uncoupled barrier. Contrary to this expectation, our calculations show that determining 
the weight factors at the barrier position does not give proper weight factors even though 
the weight factors do not have a strong radial dependence. The situation would be much 
more serious in realistic calculations, since the weights are still changing at the barrier 
position due to the fact that the coupling extends outside the Coulomb barrier JTl|. 

The same calculations were repeated for different values of the excitation energy e and 
the results are shown in Fig. 3. The quantities shown in each panel are the same as those 
in Fig. 2 except for the third panel, where only the optimum weight factors for the lower 
barrier are plotted. We again observe that the weight factor changes only marginally as a 
function of the incident energy, even when the excitation energy is finite. We thus conclude 
that the eigen-channel approach is still applicable even when the excitation energy of the 
intrinsic motion is finite. 

When the excitation energy e is zero, the barrier distribution, i.e. the first derivative 
of the penetrability has two symmetric peaks and the weight factor has no incident energy 
dependence (the solid line). As the excitation energy increases, some strength is transfered 
from the higher peak to the lower peak, and the barrier distribution becomes asymmetric. 
As the excitation energy significantly exceeds the curvature of the bare barrier, which 
is about 4 MeV in our example, one expects to reach the adiabatic limit. To illustrate 
this, figure 4 shows the influence of the coupling to an excited state whose energy is 8 
MeV. The figure also contains the result for the the no coupling case (the dotted line) 
for comparison. In this case, the weight for the higher peak is considerably smaller than 
that for the lower peak, and the barrier distribution has essentially only a single peak. 
The peak position, however, is shifted towards a lower energy. This is consistent with 
the adiabatic picture, i.e. the main effect of the coupling to a state whose excitation 
energy is much larger the barrier curvature is to introduce an energy independent shift 
of the potential accompanied with the mass renormalization |I2], [lf|. Correspondingly, 
the barrier distribution is only shifted without significantly changing its shape unless 
the coupling form factor has a strong radial dependence ||17|| . It is thus clear that the 
second derivative of Ea can represent a wide range of situations between the two extreme 
limits, i.e. from the adiabatic limit where the coupling leads to an adiabatic potential 
renormalization, to the sudden limit where the coupling gives rise to a barrier distribution. 

In summary, we have discussed the energy dependence of the weight factors in the eigen 
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channels approximation. Performing exact coupled-channels calculations, we showed that 
the energy dependence is very weak, regardless of the excitation energy of the intrinsic 
motion. Therefore, it makes sense to call the second derivative of Ea as the fusion barrier 
distribution, even when the excitation energy is finite. We also discussed a transition 
from the sudden to the adiabatic tunneling limits, and showed that the fusion barrier 
distribution can represent these two limits in a natural way. 
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APPENDIX: EIGEN-CHANNEL APPROXIMATION IN PATH 

INTEGRAL APPROACH 

Consider a system where a macroscopic degree of freedom R couples to an intrinsic 
degree of freedom £. We assume the following Hamiltonian for this system; 

h 2 d 2 

H (R, 0 = ~Y gjp + U ( R ) + H ««) + 0 , ( 8 ) 

where p is the mass for the macroscopic motion, U(R) the potential in the absence of 
the coupling, R 0 (£) the Hamiltonian for the internal motion, and V(R,£) the coupling 
between them. The barrier transmission probability from the initial position R t on the 
right side of the barrier to the final position Rf on the left side is then given by 

P(E) = lim r dT e {i,n)ET [°° dT e - {i/h)Ef 

R i-°° \ U? J J 0 Jo 

x J V[R(t)\ J P[R(t)]e (i/a)[5t(fl ’ T) - St(S ' ? VM(R(t),T;i?(t),T), (9) 

where E is the total energy of the system, and Pi and Pf are the classical momenta at 
Ri and Rf, respectively. S t (R,T ) is the action for the macroscopic motion along a path 
R(t), and is given by 

S t (R, T ) = jf dt <^-pR{tf - U{R{t))^j . (10) 

The effects of the internal degree of freedom are included in the two time influence func¬ 
tional p M , which is defined as 

p M (R(i),T-R(t),T) = J2< ni\u\R(t),T)\n f >< n f \u(R{t),T)\rii > . (11) 

n f 
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Here, n l and rif are the initial and the final states of the internal motion at position Ri 
and Rf, respectively. The time evolution operator u(R(t),t) of the intrinsic motion obeys 


ih^u(R(t),t) = [H 0 {£) + V{R{tU)]u(R{t),t). 

with the initial condition u(R,t = 0) = 1. 

We introduce the eigen-channel (or the adiabatic) basis, as 


(tfott) + V(R,0)MR,0 = K(R)MR,0- 


( 12 ) 


(13) 


Expanding the intrinsic wave function at time t in this basis 

u(R,t)\rii >= ^)a n (£)exp ~ [ dt'\ n (R(t')) 
n L n Jo J 

Eq. (2) for the barrier penetrability reads 


I <Pn(R(t)) >, 


(14) 


P(E) = lim 


ILL 

h 2 


dT e {ilK)ET / d f e -{i/h)ET 


x fv[R(t)] f V[R{i)]Y.an{T)a m {Ty < (f m (R(f))\ip n (R(T)) 

** ** n,m 

| / 0 T dt(±iJ,R 2 -U(R)-\ n (R)) ~lifo di (j» R -U{R)-X m (R)^ 


> 


x e 


(15) 


We use the semi-classical approximation, and for energies well below the barrier, where 
the single paht dominates, evaluate the path integral along the classical path 

R(t) = m = R d (t), T = T* = T cl . (16) 

In this case, the orthogonalitiy of the adiabatic basis leads to 

P(E) = E v n (E)P 0 (E- U(R) + A n (R)), (17) 

n 

where the weight factors are given by 

v n (E) = \a n (T cl )\ 2 . (18) 


The weight factors depend implicitly on the energy E through the time evolution of the 
intrinsic system along the classical path. 
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Figure Captions 


Fig. 1 : The eigen-potentials (the upper panel) and the associated weight factors (the 
lower panel) as a function of the relative distance R. The solid and the dashed lines 
in the lower panel are the weight factors for the lower and the higher eigen-potentials, 
respectively. The dotted line in the upper panel represents the bare potential barrier. 
Fig. 2 : The penetrability (the first panel) and its first derivative (the second panel) for 
the two level problem as a function of the incident energy. The third panel shows the 
optimum weight factors obtained according to eq. (5). The solid and the dashed lines 
correspond to the weight factors for the lower and the higher potentials, respectively. 
Fig. 3 : Same as fig. 2, but for different values of the excitation energy of the intrinsic 
motion which are denoted in the inset. The third panel shows only the weight factors for 
the lower potential. 

Fig. 4 : Effects of the coupling to a high excitation energy state. The dotted line is the 
result without channel coupling, while the solid line is obtained by taking into account 
the coupling to an excited state whose excitation energy is 8 MeV. 
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